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1 Introduction 



Experiments show that cellularisation of flames results in an increase of their prop- 
agation speed. In order to understand and exploit this phenomenon, we study the 
evolution of flame fronts governed by the Sivashinsky equation 

9*$ - {d^^f = 9,.,$ - (7/2)a,H[$] + /(x, t), -oo < a; < oo, t > 

(1) 

with the force term f(x,t). Here ^{x,t) is the perturbation of the plane flame 
front, 7i[$] = 7r-i£^(x - ^)^^^{^,t)d^ is the Hilbert transformation, and 
7 = 1 — ph/ pu is the contrast in densities of burnt and unbumt gases pb and 
Pu respectively. Initial perturbation 0) is given. 

The equation without the force term was obtained in ID as an asymptotic 
mathematical model of cellularisation of flames subject to the hydrodynamic flame 
instability. The force term was suggested in [2J in order to account for the effect of 
the upstream turbulence on the flame front. It is equal to the properly scaled turbu- 
lent fluctuations of the velocity field of the unbumed gas. In |l3l and H equation 
(dl) was further refined in order to include effects of the second order in 7. How- 
ever, as mentioned in [4J, this modification can be compensated upon a Galilean 
transformation combined with a nonsingular scaling. Thus, we have chosen to 
remain within the first order of accuracy in 7 of the original Sivashinsky model 
dU) as it should have the same qualitative properties as the more quantitatively 
accurate one. 

The asymptotically stable solutions to the Sivashinsky equation with f(x, t) = 
corresponding to the steadily propagating cellular flames do exist and are given 
by formula 

N 

^N,L{x,t) = VN,Lt + 2^1n I cosh27r6„/L - cos2ttx/L\, (2) 

n=l 

discovered in |5|. Here, real L > and integer N from within the range < 
< Nl = cei\{'jL/8iT + 1/2) — 1 are otherwise arbitrary parameters. Also, 
Vl = 2TtNL~^ (7 — ATtNL^^), and 61, 62, • • • , ^tv satisfy a system of nonlinear 
algebraic equations available elsewhere. Functions Q have a distinctive set of N 
complex conjugate pairs of poles Zn = ±ibn, n = 1, . . . , N and are called the 
steady coalescent pole solutions respectively. 

The steady coalescent pole solutions Q with the maximum possible number 
= Nl of the poles were found to be asymptotically, for t ^ 00, stable if the 
wavelength of the perturbations does not exceed L, see |6|. However, in spite of 
their asymptotic stability, there are perturbations of these solutions which can be 
hugely amplified over finite intervals of time resulting in significant transients, see 
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|I3. These perturbations are nonmodal, because they cannot be represented by the 
single eigenmodes of the linearised Sivashinsky equation. In what follows we are 
interested in solutions Q with = Nl and retain the index L only. Also, in all 
reported calculations 7 = 0.8. 

In this work we calculate the most amplifiable nonmodal perturbations to the 
asymptotically stable cellular solutions of the Sivashinsky equation and use them 
to investigate the response of the flame front to forcing. In particular, we study the 
effect of stochastic forcing or noise. The investigation of the effect of noise in the 
Sivashinsky equation was carried out numerically and the observations were rein- 
forced by the analytical analysis of an approximation to the linearised Sivashinsky 
equation suggested in |i8J. 



2 The largest growing perturbations 

Substituting <l>(x, t) = t) + t) into (HJ) for /(x, t) = and linearising 

it with respect to the L-periodic perturbations t), one obtains 

dtcj) = {d,^L)dA + d..(l> - {l/2)d,H[ct>\ = Al<P, 

(3) 

0(x,O) = 0) - 0). 

The operator generates the evolution operator e*^^ , which provides the solu- 
tion to Q in the form t) = e*^^0(a;, 0). 

Assuming that the polar decomposition of the evolution operator does exist, 
we write it as 

e*^^ =W(t)5(t), (4) 

* 1 /2 

where U{t) is a partially isometric and S{t) = [(e*^^)* e*^^] is the nonnega- 
tive self-adjoint operator, see e.g. f9\. The partial isometry of U{t) implies that it 
preserves the norm when mapping between the sets of values of (e*^^ ) * and e*^^ , 
i.e. ||W(t)0|| = \\4>\\. Then, under certain conditions, = 0)|| 

and for the 2-norm the sup {||5(t)0(a;, 0)|| x||0(s, 0)||"^} is equal to 

(/>(a:,0)eI'(e*'*L) 

the largest eigenvalue oi{t) of S{t). This eigenvalue is associated with the eigen- 
vector t) of S{t). 

The eigenvectors ^q(x, t) of S{t) are mutually orthogonal at any given time 

t = t* and can be used as a basis in the space of the admissible initial conditions 

00 

0) = ^ Ca{0,t*) xil>a{x,t*). Then, the associated eigenvalues o-a{t*) pro- 

a=l 

vide the magnitudes of amplification of the ipaix,t*) components of the initial 
condition (f){x, 0) by the time instance t*. Note, that for Q the 2-norm of the per- 
turbation t) is just its energy and that the eigenvalues aait), a = 1,2, .. . and 
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eigenvectors ipa{x, t) of S{t) are the singular values and the right singular vectors 
of e^^^ respectively. 

According to lITOll . the Fourier image of the operator Al is defined by the 
(fc, /)-th entry of its double infinite (—00 < k,l < 00) matrix 

{A,),,= (-^e + ^\k\)6,,+ ^-^ls^^^^^ (5) 

^ ^ n=l 

where 5k,i is the Kronecker's symbol. By limiting our consideration to the first K 
harmonics, we approximate our double infinite matrix Al with the {2K + 1) x 

{2K + 1) matrix A^^^, whose entries coincide with those of Al for —K < k,l < 

K. Then, the matrix 6*^^- ^ ~ e*^^ can be effectively evaluated by the scaling and 
squaring algorithm with a Pade approximation. Eventually, the required estima- 
tions of aa{t) and Fourier images of ^/'q,(x, t) can be obtained through the singular 

value decomposition (SVD) of 6*^^- ' , see e.g. [ 1 1 1. 

Indeed, if the SVD of e l is given by 

e'^^'^"^ = W{t)V{t)V{ty, (6) 

where W(t), V{t) are unitary and X'(t) is the nonnegative diagonal matrix, then 
the matrices 

u{t) = w{t)v{ty, s{t) = v{t)v{t)v{t)* (7) 

satisfy the adequate finite-dimensional projection of the polar decomposition Q 
and the eigenvalues (Ja{t), a = 1,2, .. . and eigenvectors 'ijjaix,t) of S{t) are 

just the singular values and the Fourier syntheses of the right singular vectors of 

-7-' (A-) 
e ^ respectively. 

Graphs showing dependence of a few largest singular values of e*"^^ versus 
time are shown in Fig. [B One may see that values of cri,2(i) for large enough 
t match the estimation of the largest possible amplification of the perturbations 
t) obtained in by a different method. An even more impressive observa- 
tion is that the dimension of the subspace of the significantly amplifiable pertur- 
bations is very low. Perturbations of only two types can be amplified by about 10^ 
times. 

The initial conditions 0), which would be the most amplified once by 
t* = 100, 200, 300 and 10^ i.e. t*), are depicted in Fig. 13 The dominating 

singular modes ipaix,t) stabilize to some limiting functions for t > 300. For 
example, their graphs for t = 500 and t = 10^ are indistinguishable in Fig. |2l 
However, they vary in time significantly when t < 300 and for t = 200 the 
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Figure 1 : Twelve largest singular values of e' 



associated amplification ai 2(200) is already about 10'^, though ^'1,2(3^, 200) does 
not coincide with neither ^^1,2(2;, 10'^) nor ^'3,4(3^) 10^). Thus, the dependence of 
ipa on time makes the dimension of the subspace of perturbations, which can be 
amplified say about lO'^ times much higher than two in contrast to what could 
be concluded from the graphs in Fig. [B This illustrates the complicatedness of 
studies of the effect of transient amplification on short time scales t < 300. 

Fourier components of t*) and ^^2(2;, t*) fort* = 10'^ are depicted in Fig. 



Evolution of the perturbations, which grow the most and is governed by the 
nonlinear Sivashinsky equation, is illustrated in Fig. IH All the profiles were 
displaced vertically in order to compensate for steady propagation of flames in 
such a way that their spatial averages are equal to zero. Matching graphs of the 
spatially averaged flame propagation speed 



are shown as well. The initial conditions were 0) = 0) + eipaix, t*), 

where e = ±10"'^, a = 1,2, and t* = 10^. The computational method used in 
this work was presented in [il2ll . 

The asymmetric singular mode ipi{x, t*) results in appearance of a small cusp 
to the left or to the right from the trough of 0) depending on the sign of 

e. After the cusp merges with the trough, the flame profile converges slowly to 
+ Ax,t), where sign(£:)Ax > 0. For a positive e = 10^^ the effect is 
illustrated in Fig. IH Graphs of t) for e = — lO^'^ are exact mirror reflections 
of those depicted in|4ta) and graphs of < $t > are exactly the same. 



m 



L/2 




(8) 
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The symmetric singular mode ip2{x, t*) produces two symmetric dents moving 
towards the trough on both sides of the profile if e < 0, see Fig. Etb). By t ^ 500 
the flame profile returns very closely to ^^(x, t). For e > two small cusps 
move towards the boundaries of the computational domain creating a quasi-steady 
structure shown in Fig. Uc) for t = 270. This structure survives until t ~ 1800, 
but eventually bifurcates, see Fig. |4td), and the solution converges to $^(0; + 
Ax,t), Ax < 0. It looks like the bifurcation in question is associated with the 
lack of the asymptotic stability of the intermediate quasi-steady structure. As 
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Figure 3: Fourier coefficients of ipa{x, t*) = i'a,k{t*)e''^''^'^^^ for t* 

k=—oo 



such, it was triggered by a random perturbation and could equally result in the 
displacement of the limiting flame front profile into the opposite direction Ax > 
0. 

Behavior of perturbations ^^1,2(2;, t*) of the amplitude e = 10^® was not as 
impressive, but they managed to produce a visible effect on the flame front profile. 
The same can be said about ^3,4 (x, t*) of the amplitude e = 10^^. Perturbations 
corresponding to V^q,(x, t*) of higher orders did not grow significantly and did not 
cause any noticeable changes to $l for e up to 10~^. 

Thus, the singular modes ^^1,2(3^, t*) should be responsible for the interaction 
of the flame front t) with all the perturbations of small enough amplitude. 
The time scale of these interactions is about 300 for L = AOn and is of order 0{L) 
in general. More singular modes ^/'q,(x, t*) of higher orders a > 2 are becoming 
important as the amplitude of the perturbations grows. The time scale of evolution 
of (j){x,t) for 0) = ipaix^t*) lessens as a grows necessitating to take into 
account the dependence of ipaix,t*) on t* and creating further problems in the 
efficient description of the subspace of important perturbations. Therefore, there 
is a critical perturbation amplitude beyond which the representation of /(x, t) in 
terms of the singular modes ^/'q (x, t* ) is not as beneficial as for smaller amplitudes . 



3 A simplified linear model 

Prior to experimenting with dU) we consider a simplified linear model suggested 
in lIHl. The L-periodic steady coalescent A^^^^-pole solution Q has a characteristic 
wavy or cellular structure and can be represented in a vicinity of the crest as 
$l(x, t) fa $L(0,t) - x^/(2_R) + 0{x^). Here, R is the radius of curvature of 
the flame front profile in the crest. For large enough L, it can be approximated as 
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Figure 4: Solutions to (HJ) for <l>(x, 0) = + ^i'ai.x, t*) and L = AOn. 



R = C1L + C2, where ci and C2 are some constants. Note, that in the approximation 
of ^L{x,t) the origin x = was chosen exactly in the crest of t). Thus, 
dx^L ~ —x/R + O(x^) in a vicinity of x = and the equation suggested in lIHl 
can be written as 

dt(j) + R-^xd^ = dxx<P + (7/2)9,7^[0] + /(x, t), (9) 

where —00 < x < 00 and t > 0. The latter equation is much simpler than Q, 
yet it is meaningful enough to study the development of perturbations of t) 
appearing in the crest. 

Equation Q can be solved exactly. Applying the Fourier transformation we 
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obtain 



dtT[<P\ - R-'^d^n<P] = - (4vr'e' - vTTiel - R-') -^[0] + nfM, t), (10) 

which is a linear non-homogeneous hyperbolic equation of the first order. Using 
the standard method of characteristics its exact solution can be written as follows 

t 

miU) = 6^(e,t)^[0^°)](le|e*/^)+ 1 g{U-r)m Oe|e(*-^)/«,r] dr, (11) 



where 

and jF[/](^,t) = f[x,t)e~'''^'^^^dx denotes the Fourier transformation of 

If the initial condition is a single harmonics cp^'^^x) = cos(27r^o2^ + ^) and 
f{x,t) = 0, then 

<jy{x, t) = e-2.^i?(i-e--/-)5g+.,ij(i-e-/H)5„ (27reoxe-*/^ + . (13) 

The infinite time limit of dTSt is equal e~^'^^^o+'^T^?o cosy) and is reached ef- 
fectively on the time scale of order 0{R). This time limit attains its maximum 
g7 ^/8cos(/9 for io = = 7/(47r), matching the asymptotic estimation of lE). 
Note that the wave number of the largest Fourier component k* of both t*) 
and ^2{x,t*) for t* > 300 is equal to ^ = k*/L = 7/(47r) as well, see Fig. 13 
A few graphical examples of function (O are given in Fig. |5I Note that the ar- 
gument of the cosine in (fT3t depends on time, which means that even if the initial 
condition (p^^^ (x) is a linear combination of mutually orthogonal cosine harmon- 
ics, then the solution (j)(x, t) will remain a linear combination of cosine harmonics 
for t > 0, but those harmonics will no longer be mutually orthogonal. This ex- 
plains why the most amplified perturbations are formed by linear combinations of 
a few initially orthogonal harmonics and approximate ^/^i 2(x, if:*) asymptotically 
for L ^ oo. 

Behaviour of ([T3t is in a sharp contrast with the evolution of the single har- 
monics perturbations of the plane flame front 

0(a;, t) = e(-^"'«o+-7€o)< cos(27r^oa; + <^), (14) 

which grow infinitely if < 7/(4vr) or decay otherwise. They are governed by 
the equation associated with a self-adjoint differential operator, which is obtained 
from ^ upon removal of the term R^^xdx4>- Solution (fT4l) does not result from 
(fT3t for R ^ oo, but is only equivalent to it when t/R <^ 1. The difference be- 
tween ([T3t and (fT4l) is an explicit illustration of the nonnormality of Q introduced 
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(a) eo = 7/(27r), = (b) = 7/(47r), if = 




(c) ^0 = 7/(87r), <^ = (d) ^0 = 7/(327r), <^ = 




(e) Co = 7/(4vr), <^ = 7r/4 (f) = 7/(4vr), = 7r/2 



Figure 5: Examples of solutions (fTSt for R = 146.7126, which corresponds to 
L = 407r. 

by the non- self adjoint term R~^xdx(l). Flattening of the crests of cellular flames 
and increasing local resemblance with the plane front as R increases was noticed 
long time ago, prompting a hypothesis of a secondary Darrieus-Landau instabil- 
ity. Model Q indicates that the hypothesis is unlikely to be correct. Although, 
because of the flattening of the crests of the flame front profile, perturbations of 
the front can be transiently amplified at a rate rapidly increasing with R, this tran- 
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sient amplification is entirely different from the infinite growth of perturbations 
in the Darrieus-Landau instability of plane flames. Moreover dynamics of pertur- 
bations in the case of cellular flames does not converge to that of the plane ones 
continuously in the limit i? — > oo. 

Solution dni), (O for 0) = e'^^ ,p> can be represented in a closed 
form as well. Routine integration yields 



pa 



nbx ^ 
cos h K 



iTTbx ^ I h + i2'Kx\ 

e « eri — — 

2^ ) 



(15) 



where 



a{t) = 27c'^R {e 



2t/R 



1) + vr^e^*/^ 



/p, b = h{t) = TT-fR (e 



t/R 



(16) 

The result is illustrated in Fig. |6l where case (c) corresponds to the maximum 
growing perturbation of type (f){x, 0) = e~^^ and initial condition 0) = 5{x) 
was used in (d). The solution formula in the latter case is given by dTSt with p for- 
mally replaced by vr and it also should be used with a = a(t) = 2u:'^R (e^*/^ — l) 
and h exactly the same as in ([T6b . 

The steady coalescent pole solutions to the Sivashinsky equation correspond 
to the flame fronts propagating steadily with the velocity exceeding mj, = 1 by 
Vl, see e.g. lfT3ll . Addition of the perturbation t) results in a change in the 
velocity of propagation by the value of the space average of dtcf), which we denote 
<(pt>- The correction provided by the t) is only valid in a small vicinity 
of the crest of t) . In sequel, the correction of the speed <0t > is only 

valid for a small region —e<x<eof the flame front in a vicinity of the crest 
of t). Hence, for our simplified linear model we define the increase of the 
flame propagation speed as follows: 



<0f>-- 



— I dtcpdx 



(17) 



For the single harmonics solution ([T3t the expression for <0t> is obvious and 
is illustrated in Fig. IT] These graphs demonstrate high sensitivity of < 0t > to 
the wavelength of the perturbation. The phase, or location of the perturbation, is 
important as well. 



4 The effect of noise 

According to the results of Section |2l the forcing in the Sivashinsky equation can 
be decomposed into the most amplifiable nonmodal component and the orthogo- 
nal complement. The latter can be neglected reducing spatio-temporal stochastic 
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Figure 6: Examples of solutions dTSt . (fT6b for i? = 146.7126, which corresponds 
to L = 407r. 
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Figure 7: Averaged increase of the local flame propagation speed <(pt> for solu- 
tions (fT3t . (y9 = 0. 



noise to the appearance of a sequence of the most growing perturbations ipa^ (x, t*). 
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1 < ttm < a* = a*(/o) at a set of time instances tm, m = 0, 1, 2, . . .: 

oo 

/(a;, t)^foYl ^""^ ^*)'^(^ - ^-)' 1 < «m < a* = a*(/o). (18) 

m=0 

Thus, the amplitude of noise /o, alongside with the averages and the standard 
deviations of tm+i — tm and am, m = 0,l,2,...are the only essential parameters 
of such a representation of noise. 

The impulse-like noise (fTSt is used here for the sake of simplicity. Some 
arguments towards its validity were suggested in |2|. More sophisticated and 
physically realistic models of temporal noise characteristics can be used with (H]) 
as well. 

If /o ^ then noise is not able to affect the flame at all and can be 

completely neglected. This case can be referred to as the noiseless regime. On the 
other hand, if /o is comparable with the amplitude a of the background solution 
t), then almost all components of noise will be able to disturb the flame 
and the f{x,t) in (HI) should be treated as a genuine spatio-temporal stochastic 
function. This is the regime of the saturated noise. 

Eventually, there is an important transitional regime when the noise amplitude 
/o is at least of order of erf ^(t*), but still much smaller than a. In this case only the 
disturbances with a significant component in the subspace spanned by the linear 
combinations of ijja{x,t*), 1 < a < a* have a potential to affect the solution. 
All other disturbances can be neglected and the force f{x, t) in the Sivashinsky 
equation dU) can be approximated by dTSt with a finite value of a*. We would 
like to stress that though such representation of noise is correct for noise of any 
amplitude, apparently it is only efficient if /o < a^}{t*) -C a, where a* is small 
enough. 

4.1 Noise in the linear model 

A random point-wise set of perturbations uniformly distributed in time and in the 
Fourier space is a suitable model for both the computational round-off errors and 
a variety of perturbations of physical origins. We are adopting such a model in 
our analysis in the following form 

Mit) 

f{x, t) = ^am cos(27r^„x + ipm)S{t - tm), (19) 

m=l 

where am, tm, ^m, and ipm are non-correlated random sequences. It is assumed 
that ti < t2 < ■ ■ ■ < tm < ■ ■ ■ < tM{t) < t, < Lfm < ^71, and ^m > 
0, m = 1, 2, ... , M{t). Availability of the exact solution STSh makes it also 
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possible to study an alternative noise model based on elementary perturbations 



\2 

^' , which are local in physical space. 



Using (fTTT) . (fT2ll for the zero initial condition, the exact solution to Q, (fT9t 
can be written as 



M(i) 



Clrn.C 



-27r2/?[l-e-2('-'™)/fl]§^+7r7_R[l-e-(*-*™)/«]5„ 



m=l 



X COS [27r^^e-(*"*'")/-^a; + (/?^] . (20) 

The expression for <(j)t> is obvious, see dTTt . and is illustrated in Fig. [H Here 
we generated random sequences of the time instances tm with a given frequency 
F = M(T) /T on a time interval t E [0, T]. Values of the wave number and of 
the amplitude Um were also randomly generated and uniformly distributed within 
certain ranges. According to the formula for <(l)t>, the effect of the phase shift 
Lfm just duplicates the am- Therefore, its value was fixed as = 0. 

If values of dm are uniformly distributed in [—1, 1], then the time average of 
<(pt> is obviously zero, because of the linearity of the problem. In the Sivashinsky 
equation this effect is compensated by the nonlinearity. The cusps generated by 
the perturbations of opposite signs move into opposite directions along the flame 
surface, see Section |2l though they both contribute into the speed positively. This 
effect of the nonlinearity can be mimicked by restricting the range of possible 
values of the amplitudes, e.g. a„i G [0, 1], as this can be seen in Fig. [U 




xlO 



xlO 



(a)F = 1, am,U e [0,1] 



(b)F = 1/33, am e [0,1], 
U = 7/(47r) 



Figure 8: The effect of noise (fT9l) on <(pt> for L = AOtt. 

Figure Etb) shows that the increase of <0f> seen in Fig. Uta) can be matched 
by using only the largest growing perturbations with much smaller frequency, 
which is quite expected in virtue of the linearity of the problem. The amplitude of 
fluctuations in Fig. [8ta) is noticeably less than in Fig. Etb). This is attributed to 
the smoothing effect of the less growing perturbations. 
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Because of the linearity of the problem in question, the effect of F = M(T)/T 
and L on <0(> is straightforward. In particular, the value of <(pt> raises up to 
about 4 X 10® for L = SOvr and other parameters the same as in Fig. Eta). It should 
be noticed however that because of the limitation > the quantity <(j)t> does 
no longer represent the increase of propagation speed of the flame, but is just a 
measure of the rate of transient amplification of perturbations. Figure |9lillustrates 
the point, see also Q and I.12J . 
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Figure 9: The effect of noise ([T9b on for L = SOvr, F = 1, G [0, 1], and 

e-e [0,1]. 

Direct studies of the effect of noise in the Sivashinsky equation necessitate use 
of numerical simulations. However, because of the intrinsic discontinuity of noise, 
such DNS are hampered with very low accuracy of approximations questioning 
the validity of numerical solutions. In this work we used explicit solutions (l20b 
in order to validate DNS of Q and, in sequel, of O- The DNS of Q, was 
carried out using a spectral method. The delta function 5{t—tm) was approximated 
by (7rAt)"^/^e"(*"*")^/^* with a small enough At ^ The calculations have 
shown that discrepancies between (l20b and its numerical counterparts obtained 
with the same sets of tm, am, and niight be noticeable in a neighbourhood of 
the time instances t tm- Although, the averaged characteristics like <(pt> were 
quite accurate. So, this linear model validates the DNS of the forced Sivashinsky 
equation at least in relation to the averaged flame propagation speed. 

4.2 The Sivashinsky equation 

We carried out a series of computations of ([l]), (fTSt with ^^(x, t) as initial con- 
dition and with a variety of parameters of the noise term. Up to twelve basis 
functions -ipaix, 10^), where a was uniformly distributed in the interval 1 < a < 
a* < 12, were used. The sign of /o in dTSt was either plus or minus for every m 
with the equal probability 1/2. The delta function 5{t — tm) was approximated by 
(7rAt)"^/^e"(*"*™)^/'^* with a small enough value of At. 
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The effect of the amplitude of noise on the flame speed is illustrated in Fig. 
[TOl Use of only two basis functions ^^1,2(2;, 10^) gives almost the same result. 
Similar to the linear model Q the only noticeable difference was in slightly larger 
fluctuations of Examples of the effect of composition of noise on 

are given in Fig. fTTltoo. 




2000 4000 6000 8000 10000 

t 



Figure 10: The effect of the amplitude of noise dTSt on <$t> for L = 407r. Here 

F = 1/15 and a* = 12. 

It was mentioned in the previous section that the wave number of the largest 
Fourier component of ^1,2(2^, 10^) is exactly the same as the wave number ^0 = 
7/(47r) of the largest growing single harmonics solution to Q. We tried to ex- 
ploit this observation and simplified (fTSt even further, replacing '4'a{m)ix,t*) by 
cos(7x/4), which corresponds to i-e. 

00 

/(x,t)^/ocos^^5(t-U. (21) 

m=0 

This kind of forcing is able to speed up the flame, but the difference between 
the computational results obtained with (fTSb and d^TTl is noticeable. It does not 
disappear even if eight nearest sidebands are added to (|2T1) . see Fig. [TTJ 
The time averages of denoted here as 

«$,»= ^ [ <<^^>dt, (22) 

tend ~ to J 
to 

are depicted in Fig. [T2lversus F (left) and /o (right). Discrepancies in for 
different a* did not exceed the variations caused by the different randomly chosen 
sequences of t^- Although, effect of using (l2Tt instead of (fTSt is appreciable. 

The correlation between the flame propagation speed and the noise amplitude 
is obvious. Note that the /o in the right most point in the graph is still about 20 
times less than the amplitude of the variation of the background solution $L(a;, t). 
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Figure 11: Examples of the effect of the composition of noise on < $t > for 
/o = 10^3 and L = 407r. 



In accordance with the idea developed in this paper, the value of <$C^t^ is 
determined by the product cti/q. It was shown in [7| that ai oc e'^^^\ resulting 
in ^$t»=«;<l>i» (e^^^^fo). Thus, the data shown in Fig. Oare at least in a 
qualitative agreement with the dependence of <s^^t^ on L, which was obtained 
in Q for a fixed noise amplitude /o ~ 10^^® associated with the computational 
round-off errors. 

Eventually, in Fig. [13] we presented the results of an attempt to control the 
flame propagation speed using our special perturbations tpa{x, t*) of properly se- 
lected amplitudes. Graphs of < $t > and <<c^t^ are shown in the left and 
numerical solution ^{x,t) corresponding to this controlling experiment is illus- 
trated in the right. The fluctuations of the obtained flame propagation speed are 
large indeed, but at least, they appear in quite a regular pattern. 

In this paper noise or forcing in ([U) represents the turbulence of the upstream 
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F fo 

(a) /o = 10-3 (b) F = 1/15, a* = 12, 



Figure 12: The effect of the composition (a) and amplitude (b) of noise on the 
spatio-temporally averaged flame propagation speed <<c;$t^. Here L = AOn and 
the temporal averaging was over the interval t E [200, 10"^]. 




Figure 13: An example of controlling the flame speed with the amplitude of the 
perturbations /q. Here F = 1/15, a* = 12, and L = AOn. 

velocity field, which is difficult to manage in practice. The controlling function 
is more effectively achieved by acoustic signals, see e.g. lfT4l . Acoustics was 
neglected in the evaluation of the Sivashinsky equation and there is no easy and 
straightforward way to incorporate it back into the model. However, because of 
a strong coupling between the velocity and pressure fields, effects of acoustic 
signals similar to those presented here can be expected as well. 

5 Conclusions 

Based on our analysis of the steadily propagating cellular flames governed by the 
Sivashinsky equation we may conclude that there are perturbations of very small 
amplitude, which can essentially affect the flame front dynamics. The subspace 



18 



formed by these special perturbations is of a very small dimension and its basis can 
be used for an efficient representation of the upstream velocity turbulence. These 
are the very perturbations which cause the increase of the flame propagation speed 
in numerical experiments. Hence, theoretically, they can be used to model certain 
regimes of flame-turbulence interaction and to control the flame propagation speed 
on purpose. 
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